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^ ; ABSTRACT 

m ■ 

• Aims. We address the issue of the multiplicity of the three brightest early-type stars of the young open cluster IC 1805, namely 
O ' HD 15570, HD 15629 and HD 15558. 

Methods. For the three stars, we measured the radial velocity by fitting Gaussian curves to line profiles in the optical domain. 
' In the case of the massive binary HD 15558, we also used a spectral disentangling method to separate the spectra of the primary 
(-H _ and of the secondary in order to derive the radial velocities of the two components. These measurements were used to compute 
^ orbital solutions for HD 15558. 

Results. For HD 15570 and HD 15629, the radial velocities do not present any significant trend attributable to a binary motion on 
time scales of a few days, nor from one year to the next. In the case of HD 15558 we obtained an improved SBl orbital solution 
with a period of about 442 days, and we report for the first time on the detection of the spectral signature of its secondary star. 
^ I We derive spectral types 05.5IH(f) and 07V for the primary and the secondary of HD 15558. We tentatively compute a first 
SB2 orbital solution although the radial velocities from the secondary star should be considered with caution. The mass ratio 
is rather high, i.e. about 3, and leads to very extreme minimum masses, in particular for the primary object. Minimum masses 
k>( , of the order of 150 ± 50 and 50 ± 15 M0 are found respectively for the primary and the secondary. 
j_j Conclusions. We propose that HD 15558 could be a triple system. This scenario could help to reconcile the very large minimum 
mass derived for the primary object with its spectral type. In addition, considering new and previously published results, we 
find that the binary frequency among O-stars in IC 1805 has a lower limit of 20 %, and that previously published values (80 %) 
are probably overestimated. 

Key words, stars: early-type - binaries: spectroscopic - stars: individual: HD 15570 - stars: individual: HD 15629 - stars: 
individual: HD 15558 
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1. Introduction 

The study of stellar populations of young open clusters 
has been the purpose of many works in the las t years (e.g. 
ISagar et al.ll988(lRaboud fc Merniilliodll998|) . One of the 

crucial questions addressed in this context concerns the 
massive star content and the binary frequency of the ear- 
liest stars harboured by these clusters. Recent numerical 
simulations of the star formation in open clusters predict 
a high level of mass segregation, along with the formation 
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of the most massive stars through simultaneous accretion 
and stellar collisions, resulting in t he presence of the mos t 
massive stars in binary systems l|Bonnell fc Batdl2002j) . 
The investigation of the multiplicity of the massive star 
population of young open cl usters is consequently of par- 
ticular interest. For instance, Inarcfa fc MermilliodI (|2nmh 
presented radial velocity measurements for 37 O- and B- 
stars in the open cluster NGC6231, and proposed binary 
frequencies for a series of open clusters, among which is 
IC 1805. 

IC 1805 is a young rich open cluster in the core of the 
CasOB6 association, in the molecu lar cloud W4 in the 
Perseus spiral arm of our Galaxy. iMassev et alJ (^2^ 
inferred an age of 1-3 Myr for the cluster, in agreement 
with other previous estimates (see iFeinstein et al.l Il986l 
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and references therein, although these latter authors 
proposed an ag e < 1 Myr) . The spectroscopic study of 
IShi fc Hi J 1119991 ) revealed that about 80 of the members of 
IC 1 805 are O or B stars. ICarcfa fc Mermill iod (2001) (see 
also Ilshidalll97(ll) estimated a binary frequency of 80% 
among the 10 0- stars in IC 1805. However, in a previous 
paper ( R,auw fc Dc Bcck^ 120041 Paper I) , we showed 
that two suspected binaries, BD-|-60° 501 (07V((f))) and 
BD+60°513 (07.5V((f))), were most probably single 
stars, suggesting that the b inary frequency proposed by 
ICarcia fc Mermilliodl l)200l|) might be overestimated. In 
this paper, we investigate the multiplicity of the three 
earhest 0-type stars of the cluster: HD 15570, HD 15629, 
and HD 15558. 



HD 15570 (04If+, V = 8.10) was proposed to be the 
most massive member of IC 1805, and incidentally one 
of the most massive and most luminous stars known 
in our Galaxy, with a present evolutionary mass of 
about 80 Mq following lllerrero et al. (2000). HD 15629 
(05V((f)), V = 8.42) has also been proposed to be a 
very massive star. According to lHerrero et all l)2000|) . this 
star could have evolved from an initial mass of about 
70 M0 , and presently be on the way towards evolutionary 
stages close to those of HD 14947 (05If+) and HD 210839 
(06I(n)fp, ACep). The same authors inferred a present- 
day evolutionary mass of about 61 Mq. For both stars, no 
clear evidence of binarity has been reported in the liter- 
ature, even though some radial velocities quoted in the 
WEBDA data base^ suggest the occurrence of variations. 

HD 15558 (my - = 8.0) was classified as an 05111(f) star 
bv lMathvsl l|1989(l . It was reported to be a spectroscopic 
binary with a period of a bout 420 d for the first time 
by Trumpler (according to lUnderhill ll967t) . Up to now, 
the onl y orbital solution p r oposed for this system is 
that of ICarmanv fc Mass'evI l)l98l|) . revealing a highly 
eccentric binary (e = 0.54) with a period of about 440 d, 
i.e. the longest of any 0-type spectroscopic binary 
known at th at time. This star i s also believed to be 
very massive. iHerrero et all (|^o3) derived spectroscopic 
and evolutionary masses of about 90 M0 for HD 15558. 
Finally, this massive binary is classified as a non-thermal 
radio emitter. In this co ntext, the study of its binarity 
is of crucial interest (see iDe BeckeilEo05l for a detailed 
discussion) . 

This paper is organized as follows. SectiorJSl presents 
the optical spectrum of the stars. The results of the ra- 
dial velocity study for HD 15570 and HD 15629 are given 
in Sect.121 The case of HD 15558 is discussed in detail in 
Sect. 01 Section|Sl consists of a discussion of our results. 
A summary of the main results and the conclusions are 
provided in Sect. El 
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2. The optical spectrum 

2.1. HD 15570 

A detailed description of our data is given in AppendixIXI 
The most prominent features observed in the spectrum 
presented in Fig. are the He 11 A 4686 and Ha lines, in 
strong emission. The He 11 A 4686 emission is one of the 
characteristics of the Of"*" type, along with Si iv AA 4088- 
4115 also in emission and Niii AA 4634-4641. We find also 
weaker emission lines of Niv at 4058 A, Cm at 4650 and 
5696 A, and also Civ A 4662. The weak emission lines 
at 6721, 6728 and 6731 A might be attributed to Cm 
llWalbornll2(¥)l . We report also on the probable presence 
of Si IV A 4656 and N v AA 6716,6718 in emission l|Walbornl 
I2OOII) . The H/3 line appears as a P Cygni profile with a 
very deep absorption component. Other absorption lines 
are mainly due to He II. He i A 4471 is very weak. We also 
note that the He i A 5875 line appears as a P Cygni profile, 
even though the emission component is rather weak. The 
Mgii A 4481 line is detected, and we also note the pres- 
ence of the Nv AA 4604,4620 lines in absorption. Even 
though its presence is not obvious in the mean SPM spec- 
trum (see Fig.nj, a weak emission feature at about 4487 
A is observed in OHP spectra and is probably due to S iv 
l|Werner fc Rauch.200H) . 

Us ing the classification criterion given by iMathvsl 
l(l988l) . we infer the spectral type 04. The various emis- 
sion lines discussed above lead to the 04If+ spectral 
type for HD 15570, in agreement with Walborn's classifica- 
tion CWal borniliiT^ . A co mparison with the spectrum o f 
the 04If+ star HD 269698 jWalborn fc Fitzpatrickl[l990|) 
lends further support to this classification. 



2.2. HD 15629 

The strongest lines in the blue spectrum of HD 15629 
(Fig. 121) are in order of decreasing intensity: H/3, Hen A 
4542, He II A 4686 and Hei A 4471, all in absorption. The 
He I A 4713 absorption line is very weak. The Mgii A 4481 
absorption line is also present. The N ill AA 4634-4641 lines 
are the strongest emission features observed in this spec- 
trum. The C III blend at about 4650 is also present and we 
note the possible presence of C iv A 4662 in weak emission. 
The weak e mission feature at abou t 4487 A is probably 
due to Siv (| Werner fc Rauchll200l|) . 

Adopting the classification criterion porposed by 

iMathv s ( 1988), we derive an 05 spectral type. Considering 
that the Nili AA 4634-4641 features are in emission, and 
that Hell A 4686 is in strong absorption, we adopt the 
05V((f)) spectral type for HD 15629, in agreement with 
previous classifications. 



2.3. HD 15558 

In Fig. 13 we present the spectrum of HD 15558 between 
3910 and 6750 A obtained with the Elodie spectrograph 
on 7 January 2004, i.e. close to the maximum of the radial 
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Fig. 1. Mean spectrum of HD 15570 calculated from spectra obtained at SPM between HJD 2 453 287.93 and HJD 
2 453 290.00. 
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Fig. 2. Mean normalized spectrum of HD 15629 between 4465 A and 4890 A as observed in September 2002 at OHP. 



velocity curve presented in Fig.01 We note that the orders 
of the echelle spectra obtained with Elodie are rather nar- 
row, rendering the rectification procedure somewhat dif- 
ficult where the wings of several lines merge over several 
tens of A . So, for the purpose of the presentation of the 
optical spectrum of HD 15558, the spectral domain be- 
tween 4580 and 4710 A in Fig.|21 comes from one of the 
Aurelie spectra obtained during the September 2001 ob- 
serving run, i.e. very close to the orbital phase of the se- 
lected Elodie spectrum. 

The strongest absorption features are the hydrogen 
Balmer lines. The Heii lines appear also in absorption. 



except for Hell A 4686 which displays a P-Cygni profile. 
The He I absorption lines at 4471 and 5875 A are also 
present, along with Hei A 4713 which is very weak. The 
other absorption features observed in the spectrum are 
Mgii A 4481, Nv AA 4604,4620, Sim A 4627, Oiii A 5592 
and Civ AA 5801,5812. The strongest emission lines are 
the Niii AA 4634,4641 lines. We report also on notice- 
able C III lines at 4650, 4652, 5696, 6721,6728 and 6731 A . 
The thr ee latter l i nes w ere qualified as selective emission 
lines bv lWalbornI l|200l|) as well as the Nv AA 6716,6718 
emission features. We report also on the presence of the 
Si IV A 4116 emission line, along with that of the much 
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Fig. 3. Spectrum of HD 15558 obtained about at phase 0.8, i.e. at the maximum of the radial velocity curve presented 
in Fig.^ This spectrum was obtained with the Elodie echelle spectrograph on 7 January 2004, except for the spectral 
domain between 4580 and 4710 A which is taken from one of the spectra obtained in September 2001 with the Aurelie 
spectrograph. We note that a large portion of the spectrum was skipped because of a lack of important spectral 
features. 



weaker Siiv A 4088. Finally, the Civ A 4662 emission line 
is cleary present. The spectral classification of the stars in 
this binary system is detailed in Sect. 15. ll 

3. The probably single stars HD 15570 and 
HD 15629 

Both stars have been mainly observed at the Observatoire 
de Haute-Provence (OHP, France) with the 1.52 m tele- 
scope. Some spectra were also obtained with the 2.1m 
telescope at the Observatorio Astronomico Nacional of 
San Pedro Martir (SPM, Mexico). For details on the obser- 
vations and the data reduction procedure, see AppendixIXI 
We have measured the radial velocity (RV) of the 
strongest lines for HD 15570 and HD 15629 by fitting 
Gaussians to the line profiles. The RVs of various lines 
from OHP data only are quoted in Table lRTI 

In the case of HD 15570, the dispersion of the RV 
of the He I A 4471 line is larger than for other lines in 
most of our data sets, as its profile deviates significantly 
from a Gaussian. In the case of the Heii A 4542 line, 
we performed a Fourier analysis of the R Vs following 
th e method des c ribed by iHeck et all l|l985fl and revised 
by iGosset et al.l ( 2001 ) . The periodogram reveals several 



peaks. The highest one is found at a frequency of 0.71 d~^. 
However, the corresponding semi-amplitude is very low 
(about 4kms^^) compared to the typical error on the 
RVs. This error, i.e. about lOkms"^, corresponds to the 
standard deviation determined for the radial velocity of a 
Diffuse Interstellar Band (DIB) at about 4762 A . The ra- 
dial velocities of the N ill emission lines appear stable on 
the time scales investigated in this study. As a result, we 
do not detect any significant RV variation attributable to 
binarity on time scales of a few days, nor from one year to 
the next, in our OHP data. The radial velocities measured 
on the main absorption and emission lines observed in the 
SPM spectra did not reveal any significant RV change nei- 
ther on a time scale of a few days, although we may note 
that the dispersion obtained in most cases is larger due to 
non-Gaussian profiles and to a poorer quality of the data. 
The lack of significant RV variations is not in agreement 
with the RV changes suggested by the data collected in 
the WEBDA data base. We caution however that these 
archive RV measurements were obtained on the basis of 
sometimes heterogeneous line lists and with instruments 
of different capabilities, at very different epochs. 

The strong profile vari ability displayed by the Heii 
A 4686 and H/3 lines (see iDe Beckeill2005|) . along with 
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their noticeable asymmetry, prevented us from obtaining 
accurate and reliable RVs from the entire profiles of these 
lines. Still, we mention that the RVs measured at the top 
(resp. bottom) of the Hen A 4686 (resp. H/3) line reveal 
rather strong wavelength shifts, which are correlated 
both in direction and amplitude (semi-amplitude of 
~30kms^^) for the two lines. However, the lack of 
significant shift for the RVs of He i A 4471, Hen A 4542 
and N ill AA 4634-4641 does not support the idea that this 
could be due to the orbital motion in a binary system. We 
finally note that the RVs of the Balmer lines measured 
on SPM spectra display the expected progression due to 
the transition from purely absorption (H9) to P-Cygni 
(H/3) profiles. 

For HD 15629, we measured the RV for several lines 
that are free from profile variability. Therefore, only the 
He I A 4471, Hen A 4542, and Niii AA 4634-4641 lines 
are considered in this discussion. A Fourier analysis of 
the RVs from the three data sets reveals a highest peak 
at 0.08 d-i (P = 12.6d) for Hei A 4471, and at 0.16 d"! 
(P = 6.25 d) for He ii A 4542. The periodogram for this lat- 
ter line presents also a peak close to that of the He i line 
at 0.09d~^ (P = ll.ld). However, the amplitude of these 
peaks is only about 5kms~^. We detect no significant RV 
variations for any of the lines. Thi s result is i n con trast 
with the variable status reported bv lUnderhilJ l(l967l) and 
[Humphrcvs (197^. We also mention that long term RV 
variations were reported for this star (see the WEBDA 
data base) , but we should consider these values with cau- 
tion for the same reasons as for HD 15570 (see above). In 
summary, we did not find any trend pointing to a binary 
scenario on the time scales sampled by our data. 

4. The massive binary HD 15558 

4.1. Radial velocity time series 

As mentioned in Sect.^ the only o rbital solution avail- 
able f or HD 15558 was proposed by ICarmanv fc MassevI 
These authors reported on a period of 439.3 ± 
1.0 d, with an eccentricity of 0.54 ± 0.05. However, their 
orbital solution was based on time series including a some- 
what heterogeneous set of radial velocities. They indeed 
used mean RVs calculated on the basis of 4 to 12 different 
lines. Among these lines, many are expected to be at least 
partly produced in the stellar wind, and might therefore 
not reflect the RV of the star itself. The heterogeneity of 
their RV time series could have a significant impact on the 
orbital parameters. In this section, we describe the proce- 
dure we adopted to establish our RV time series, and how 
we use it to determine the SBl orbital parameters of the 
system. 

The details on our data on HD 15558 are given in 
AppendixIXI First of all, we focused on the wavelength 
range covered by all our spectra, whatever the instrumen- 
tation used, i.e. the spectral domain between 4455 and 
4680 A . We then selected the lines whose profile did not 



deviate too strongly from a typical Gaussian shape. Only 
the He II A 4542 and the Niii AA 4634,4641 lines meet our 
criteria. We then determined the RVs by fitting Gaussians 
to the profiles. In the case of Elodie data, the Hen line is 
simultaneously observed in two adjacent orders. We there- 
fore measured the RV on each order and we used the 
mean of the two values in our time series. The RVs ob- 
tained from these three lines, along with those obtained 
from other lines observed only in our echelle spectra, are 
compiled in Table l5?2l In this Table, the column labeled 
'Mean' contains the mean RVs obtained from the three 
lines discussed above. We estimate that the expected un- 
certainty on the RVs are respectively of the order of 15- 
20, 10-15, and 5-8kms~^ respectively for low resolution 
Aurelie, medium resolution Aurelie, and high resolution 
Elodie data. However, we mention that the error on the 
RVs measured on some low quality Elodie spectra may be 
significantly larger. These uncertainties were estimated on 
the basis of radial velocity measurements performed on 
DIBs. 

4.2. Period determination 

We performed a Fourier analysis o n our RV tin i e serie s fol- 
lowing the technique des c ribed by iHeck et all l|l985l) and 
revised by iGosset et al.l (l200ll). as us ed for instance in 
Paper I and bv lDe Becker et al.l l)2004(l for other spectro- 
scopic binaries. We independently applied the same tech- 
nique to three RV time series obtained respectively from 
Hen A 4542, a mean of the RVs from the two Niii AA 
4634,4641 lines, and a mean of the Hen and Niii lines. 
In the three cases, the periodogram is dominated by a 
strong peak at a frequency of 0.00226 d~^, corresponding 
to a period of about 442 d. As the time base covered by 
our data is about 1608 d, the typical width of the peaks 
of the periodograms is about 6.22 x lO"'' d~^. Considering 
that the uncertainty on the frequency is about 10 % of the 
width of the peak, we obtain an uncertainty on the period 
of about 12 d. The fact that the three RV time series lead 
exactly to the same period is a strong argument for this 
value being the true period of the system. 

4.3. Orbital solution 

We obtained the SBl orbital solutio n of HD 15558 us- 
ing the method of IWolfe et alJ l)l967j) for SBl systems. 
We assigned different weights to take into account the 
expected uncertainties affecting our RV measurements. 
These weights vary between 0.1 and 1.0 respectively for 
very poor quality and good quality Elodie data. Aurelie 
RVs get intermediate weights to take into account the fact 
that though the resolution of the spectrograph is rather 
low, the quality of the data is much better than for Elodie 
spectra. 

We first fixed the period to the value determined from 
the Fourier analysis of our RV time series, i.e. ~ 442 d. We 
obtained similar results for the orbital parameters, i.e. for 
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Fig. 4. Radial velocity curve of HD 15558 for an orbital 
period of 442 d. The open and filled hexagons stand re- 
spectively for the primary RVs obtained with Elodie and 
Aurelie spectra. The solid line yields our best fit orbital 
solution respectively corresponding to the parameters pro- 
vided by Tabled The mean of the RVs of the He ii A 4542 
and of the Niii AA 4634,4641 lines was used to compute 
this orbital solution. 

Table 1. Orbital parameters of the SBl solution of 
HD 15558. To refers to the time of periastron passage. 7, 
K, and a sin i denote respectively the systemic velocity, 
the amplitude of the radial velocity curve, and the pro- 
jected separation between the centre of the star and the 
centre of mass of the binary system. The last row pro- 
vides the mass function. The orbital elements are given 
for the solutions obtained respectively from the mean of 
the RVs of the Hen A 4542 and Niii AA 4634,4641 lines, 
and from the RVs of He 11 A 4542 only. To is expressed in 
HJD-2 450 000. 





Hen A 4542 
Niii AA 4634,4641 


He II A 4542 


P (days) 


442 (fixed) 


e 


0.39 ± 0.03 


0.40 ± 0.03 


To 


1795.400 ± 6.556 


1798.385 ± 7.207 


7 (kms""*") 


-50.2 ± 1.1 


-40.1 ± 1.1 


K (kms~^) 


43.9 ± 1.8 


41.3 ± 1.8 


UJ 


116° ± 6° 


121° ± 7° 


a sin i (R©) 


350.0 ± 15.2 


331.2 ± 15.2 


/(m) (Mq) 


3.0 ± 0.4 


2.5 ± 0.3 



most of them within the error bars, whatever the RV time 
series used. We also calculated the orbital solution through 
an iterative process allowing the period to vary, but it did 
not improve significantly the results. We indeed obtained a 
period of about 445 d, which can not be distinguished from 
the 442 d period obtained in Sect l4.2l provided the uncer- 
tainty on the period is about 12 d. Therefore, we adopted 
the results obtained with a period fixed to 442 d. The cor- 
responding orbital parameters are quoted in Tabled for 
the series of RVs obtained respectively on the basis of the 



Hen and Nin lines {left column), and on the basis of the 
He II line only right column. 

4.4. Searching for the companion 

We inspected more carefully spectra obtained at phases 
close to the extrema of the radial velocity curve presented 
in Fig.d Though we did not obtain any echelle spectra 
close to the minimum of the primary radial velocity curve, 
the September 2000 and September 2001 Aurelie observ- 
ing runs fall respectively very close to the minimum and 
maximum of the radial velocity curve. We were therefore 
able to perform a more careful inspection of the line pro- 
files on the basis of our Aurelie spectra. We added together 
the 12 spectra obtained during each observing run to ob- 
tain two high signal-to-noise ratio spectra, respectively of 
about 1000 at the minimum and 900 at the maximum. We 
detected opposite asymmetries at both extrema for the 
profiles of the Hei A 4471, Mgii A 4481 and Hen A 4542 
lines, suggesting clearly the presence of the secondary. The 
inspection of the Elodie spectrum closest to the maximum 
of the primary RV curve reveals also the signature of the 
secondary in the profile of the Civ A 5812 and He i A 5876 
lines. 

The upper and lower panels of Fig.dshow the profile 
of the He I A 4471, Mgii A 4481 and Hen A 4542 lines 
respectively at the minimum and at the maximum of the 
primary radial velocity curve. We disentangled the profiles 
from the primary and the secondary by fitting Gaussians 
following an iterative process. We constructed fitting func- 
tions of the form given by Eq.d 

~^ CTj \/ ZTT L Z O j -I 

with N = 2,3,4 following the complexity of the blended 
profile. The three parameters to be determined for each 
Gaussian are respectively the normalization factor (Aj), 
the standard deviation of the Gaussian (aj ) , and the cen- 
tral position (Acj). We iteratively changed the values for 
the three parameters with N = 2 to obtain a first order fit 
of the primary and of the secondary for the helium lines. 
We then increased the number of Gaussian components in 
our fit to account for the presence of the N in A 4537 line 
on the blue side of the Hen A 4542 line. We optimized the 
fit of the He 11 line by fitting the N ill line in the blend at 
the minimum of the radial velocity curve of the primary, 
i.e. when the secondary Hen line is shifted to the red and 
consequently well separated from the Nlll line. We then 
fixed the width and the normalization factor of the N III 
line to the same values to fit the helium components of the 
blended profile at the maximum of the RV curve. We used 
a N = 4 function to fit the blend of the He i A 4471 line 
with that of Mgil at 4481 A (primary and secondary). The 
synthetic profiles obtained on the basis of Eq.d are over- 
plotted on the observed profiles in Figd While performing 
these fits, we first fixed the width and the normalization 
coefficient of the primary component to the same values 
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at both extrema. We then improved the fits by allowing 
these parameters to vary. 
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Fig. 5. Disentangling of the line profiles from the two compo- 
nents of HD 15558 respectively at the minimum of the radial 
velocity curve {upper part of the figure), and at the maximum 
of the radial velocity curve {lower part of the figure). The two 
data sets correspond to the mean spectra obtained respectively 
during the September 2000 and September 2001 observing runs 
with the 1.52 m OHP telescope. The profiles of the He I A 4471, 
the Mgll A 4481, and the Hell A 4542 lines are displayed, along 
with the best-fit functions (dotted-dashed lines). The individ- 
ual Gaussians are also displayed, with the dashed and dot- 
ted lines standing respectively for the primary (P) and the 
secondary (S). The bottom part of each panel represents the 
residuals (in the sense data - fit) . 



On the basis of our best fits, we estimated the equiv- 
alenth width (EW) of the He i and He ii lines to obtain a 
spectral classification of the two components of the sys- 
tem. At the maximum of the RV curve, we obtain EWs 
of 0.18 ± 0.02 and 0.53 ± 0.03 A respectively for the Hel 
A 4471 and He ii A 4542 lines of the primary. At the same 
phase but for the secondary, we obtain respectively EWs 
of 0.05 ± 0.02 and 0. 08 ± 0.0 2 A. U sing the classification 
criteria proposed bv iMathv^ l|l988fl . we obtain 05 and 
06.5 spectral types respectively for the primary and the 
secondary. The best-fit parameters obtained at the phase 
corresponding to the minimum of the SBl RV curve lead 
to EWs of 0.19 ± 0.02 and 0.52 ± 0.03 A (respectively 
0.08 ± 0.02 and 0.07 ± 0.03 A ) for the He i and He ii lines 
of the primary (resp. secondary). In this case, the spectral 
types of the two components of HD 15558 are 05.5 and 
07.5. Considering mean values of the EWs obtained at 
the extreme phases of the orbit, the classification of the 
two components of the system are 05.5 and 07 respec- 
tively for the primary and the secondary. 



We note that the EW of He i A 4471 undergoes a signif- 
icant decrease (~ 40%) as the secondary is receding from 
the observer. Such a decrease is also marginally observed 
for the Mg ii line, but no such behaviour is detected for the 
He II A 4542 line. However, the EWs of the primary remain 
steady from one extremum to the other. This behaviour 
is quite reminiscent of the so-called Struve-Sahade effect 
obs erved in the case of se veral massive binary system (see 
e.g. lBaenuolo et alJl999j) . However, we note that, for such 
weak lines as in HD 15558, it might be due to errors in the 
normalization of the spectra. 

4.5. Estimation of the masses of the components of 
HD 15558 

As we were able to separate the primary and secondary 
components from the profiles of a few lines at both ex- 
trema of the RV curve, we obtained a first estimate of 
the amplitude of the secondary RV curve. In the case of 
the Hell A 4542 line, we estimated that the RV of the 
secondary was determined with an uncertainty of about 
15 km s~^ at its maximum, and about 20 km at its min- 
imum. This difference is due to the fact that the blend with 
the N III line makes the determination of the parameters 
of the secondary component of the Heii line more diffi- 
cult (see bottom right panel of Fig.[31). According to our 
fit, and using the K value derived from our SBl fit (see 
TableQ]), we obtain a mass ratio of 3.3 ± 0.4. The explicit 
expression of the mass function as a function of the period 
(P), the eccentricity (e), the primary RV semi-amplitude 
{K) and the mass ratio (q) allowed us to estimate the 
minimum masses of the two components. We obtain min- 
imum masses of 152 ± 51 and 46 ± 11 Mq respectively 
for the primary and the secondary. We note that the large 
errors on the minimum masses are clearly dominated by 
the large uncertainty on the mass ratio. 

The very large masses of the two stars, mainly of the 
primary, are compatible with the large value s already men- 
tione d in the literature for this star (see e.g. lHerrero et alJ 
I2OOOI) . This is the first time that spectroscopic data lend 
support to this assertion. However, we emphasize that the 
minimum masses are poorly constrained. The critical is- 
sue regarding our results is the estimate of the mass ratio. 
More RV measurements are needed close to the extrema 
of the radial velocity curve to reduce the uncertainty on 
q, and accordingly reduce the error bars on the minimum 
masses. In order to address this issue and to tentatively 
derive the first SB2 orbital solution for this system, we ap- 
plied a distcntangling method to our spectral time series. 

4.6. Disentangling of the primary and secondary 
spectra of HD 15558 

We used the spectral dise ntangling method described by 
lOonzalez fc Levatol ( 2006|) . It consists of an iterative pro- 
cedure allowing to compute the spectra and RVs of the two 
components of a binary system. At the starting point, a 
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flat spectrum is adopted for the secondary, and the phase- 
shifted mean of the observed spectra is used as the start- 
ing primary spectrum. At each iteration, the spectra of 
the primary and the secondary obtained at the previous 
step are subtracted from the observed spectra, and the 
radial velocities are determined on the residuals through 
a cross-correlation procedure using a template including a 
series of lines. The advantage of this procedure is the fact 
that it does not require any accurate a priori knowledge 
of the spectrum of the two stars. 

As our spectral time series include data obtained with 
different instrumentations, we once again limited our in- 
vestigation to a wavelength domain covered by all our 
data. We applied the disentangling procedure to our 70 
spectra between 4456 and 4567 A, in order to include 
the He I A 4471, the Mgii A 4482 and the Heii A 4542 
lines. These lines were selected because they display rather 
clearly the signature of both components of the system. 
Starting from the two sets of RVs used to obtain the SBl 
solutions described in Tabled we applied this technique 
in three different situations: 

(a) First, we used all our spectra and we allowed both pri- 
mary and secondary RVs to vary. The RVs obtained 
after 20 iterations were used to compute an SB2 or- 
bital solution, fo llowing the technique described by 
ISana et "all l)2003(l . The orbital parameters are given 
in Table|21 The most striking results are the huge min- 
imum masses for the stars, mostly for the primary, in 
agreement with the values discussed in Sect . 14. 51 We 
see that most parameters have values close to those 
obtained for the SBl orbital solution (see TableCJ. 
Even the difference in the semi-amplitude of the pri- 
mary radial velocity curve is rather small, whilst the 
RVs were estimated on the basis of two very different 
approaches. The Gaussian fit used for the SBl case 
likely underestimates the amplitude of the RV curve, 
but Fig.El suggests that the shift between the blended 
profile and the primary component centroids should 
be very small. In the case of the mean Hell and Niii 
(resp. Hell alone) RVs, we note that the disentangling 
of four (resp. one) spectra gave unsatisfactory results 
as the derived RVs deviated strongly from those of 
other spectra obtained only a few days before or af- 
ter. The disentangling program yields indeed unrea- 
sonable radial velocities for the secondary, mostly in 
the case of spectra with poor signal-to-noise ratios (all 
deviant points come from Elodie data). For this rea- 
son, we discarded these spectra from the computation 
of the orbital solutions presented in TableEl In Fig.El 
we present the radial velocity curves obtained for both 
components of the system. We see that the agreement 
between the data and the computed RV curve for the 
secondary is rather poor, mostly at phases where the 
separation between the two components is small. The 
separated spectra of the primary and of the secondary 
in the wavelength domain used for the disentangling 
procedure are individually displayed in Fig.[7| 




Phase 

Fig. 6. Radial velocity curves of the two components of 
HD 15558 for an orbital period of ~ 442 d. The hexagons 
(resp. triangles) stand for the primary (resp. secondary) 
RVs. The solid and dashed line yield our best fit orbital 
solution respectively for the primary and the secondary, 
with the parameters provided in the left part of Table|21 
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Fig. 7. Separated spectra of the primary (lower spectrum) 
and of the secondary [upper spectrum) of HD 15558 in the 
wavelength domain used for the disentangling procedure. 

(b) We then applied the same procedure, but we fixed the 
radial velocities of the primary to the values taken for 
the SBl solution. We obtained rather similar results as 
in case (a) although the RVs derived for the secondary 
in the case of several spectra were not acceptable. 

(c) We finally applied the spectral disentangling method 
to a reduced spectral series, including only 32 spectra 
obtained at phases close to the extrema of the radial 
velocity curve, i.e. between phases 0.05 and 0.25 and 
between phases 0.75 and 0.95. We selected these spec- 
tra in order to obtain an estimate of the mass ratio un- 
affected by the data obtained when the primary and 
secondary components are heavily blended. This ap- 
proach leads to a mass ratio of 3.7 ± 1.2, but we note 
that this value should be considered with caution as 
the eccentricity (0.1 ± 0.1) and the To deviate signif- 
icantly from the values derived from the SBl and the 
SB2 (cases a and b) solutions. 

In summary, we find that the mass ratio of the two 
components of the system obtained from radial velocities 
determined on the basis of the disentangling method is 
of the order of 3, whatever the initial RVs considered as 



M. De Becker et al.: Early- type stars in the young open cluster IC 1805 



9 



Table 2. Orbital parameters of the SB2 solu tion of HD 15558 determined from the RVs computed with the disentan- 
gling procedure of lGonzalez fc Levatol l)2006|) . We started the iteration procedure with the mean of the Heii A 4542 
and the Niii AA 4634-4641 lines {left part) and of the Hell A 4542 line alone {right part). The parameters have the 
same meaning as in Tabled _Rrl sta nds for the radi us of a sphere with a volume equal to that of the Roche lobe 
computed according to the formula of lEggletonI l)l983l) . 





Hell & Niii 




He II 




Primary Secondary 


Primary 


Secondary 


P (days) 


~442 (fi: 


xed) 




e 


0.41 ± 0.06 


0.37 ± 0.07 


To (HJD-2 450 000) 


1790.377 ± 11.935 


1796.271 ± 12.970 




100° ± 10° 


120° ± 12° 


7 (kms~^) 


-54.1 ± 8.2 0.6 ± 11.4 


-51.5 ± 12.3 


-10.2 ± 15.3 


K (kms-i) 


48.3 ± 4.6 157.9 ± 15.2 


43.3 ± 6.3 


135.9 ± 19.7 


a sin i (R©) 


385.6 ± 38.8 1259.5 ± 126.6 


351.4 ± 52.0 


1101.3 ± 163.0 


m sin^ i (M©) 


234.0 ± 61.4 71.6 ± 15.8 


159.5 ± 58.4 


50.9 ± 14.4 


q = mi/m2 


3.27 ± 0.45 


3.13 ± 0.71 


iiRL sin i (Rq) 


186.5 ± 19.3 355.8 ± 37.9 


168.6 ± 26.0 


314.6 ± 52.3 



the starting point of the procedure. This value is in agree- 
ment with the mass ratio derived from the simultaneous 
fit of Gaussians discussed in Sect. 14. 51 but we failed to re- 
duce the uncertainty on q. Consequently, the minimum 
masses estimated with both approaches are similar. The 
extreme minimum masses seem therefore to be robust. 
For initial RVs obtained with the Heii A 4542 line, both 
methods provide minimum masses of the order of 150 ± 50 
and 50 ± 15 Mq respectively for the primary and the sec- 
ondary. For RVs obtained from the He ii and N ill lines, the 
minimum masses are somewhat larger with a larger uncer- 
tainty. Although the uncertainty on the secondary radial 
velocities is large, the fact that our results rely partly on 
several spectra obtained close to the extrema of the ra- 
dial velocity curves lends some support to the mass ratios 
derived by our analysis. 



5. Discussion 

5.1. Detailed spectral classification of HD 15558 

As shown bv IPe Beckeil l|2005l) , the He ii A 4686 line shows 
a P-Cygni profile for the primary and is in absorption for 
the secondary. These profiles, along with the luminosity 
contrast between the primary and the secondary, are com- 
patible with a difference of one luminosity class between 
the two stars. Considering the fact that Ha is in absorp- 
tion (see Fig.EJ, it is unlikely that the primary is a super- 
giant and we therefore propose giant and main-sequence 
luminosity classes respectively for the primary and the 
secondary. On the other hand, it is not clear whether the 
N III AA 4634,4641 lines are in emission for the secondary 
as well. Consequently, we propose that the massive binary 
HD 15558 is formed of an 05.5111(f) primary with an 07V 
secondary. 



5.2. A very massive primary? 

The most interesting - and also puzzling - result pre- 
sented in Sect. ^ is the very high minimum mass derived 
for the primary of HD 15558, especially considering its 
spectral type. On the one hand, the most massive star 
whose mass was determined using spectroscopic and pho- 
tometric data is the massive binary WR20a, with masses 
of 82 and 83 Mp, respectively for the primary and the sec- 
ondary (Rauw et al. 2004; Bonanos et al. ,2004.) . On the 
other hand, very high stellar masses were proposed (i) for 
the Pistol Star, close t o the Gal actic center, with a mass 

around 200 - 250 M© ijFieer et al 1998) , and (ii) for the 

most massive stars in the Large Mag ellanic Cloud with 
masses of the order of 120 - 200 M© (jMassev fc Hunted 
il99&) . 

The fact that we derive such a high mass for the 
primary raises the question of the so-called stellar 'upper 
mass limit'. Massive stars are expected to be vibrationally 
unstable above a given critical mass which depends on 
the evolutionary state. For zero-age main sequence stars, 
the most accurate calculations resulted in cr itical masses 
betw een about 90 an d 440 solar mas s es (se e lAppenzelleil 
Il987l for a review). lOev fc Clark3 l)2005l) statistically 
demonstrated that the local census of massive stars 
oberved so far (Milky Way + Magellanic Clouds) exhibits 
a 'universal' upper mass cutoff around 120 - 200 M© 
for a Salpeter initial mass function (IMF). Considering 
that the largest stellar mass observed in IC 1805 is 
about 100 M0, the same authors estimate that the 
probability that the stellar population of IC 1805 extends 
to 150 M© (resp. 200 M©) is about 0.63 (resp. 0.51). As 
a consequence, a stellar mass of 150±50M©, assuming 
the minimum masses derived above are close to the 
true masses {i ~ 90°), does not disagree with both the 
theoretical and the statistical approaches of the stellar 
upper mass limit. 
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Provided the mass derived for the primary is indeed 
of the order of 150 Mq^, the main issue comes from its 
spectral type. We may indeed expect the most massive 
stars to be also the earliest ones. Considering; f or instance 
the typical masses given bv lHowarth fc Prinial |l989() . we 
see that the expected masses for 05.5 stars should lie be- 
tween 50 and 80 Mq depending on the luminosity class. An 
05.5III star with a mass of about 150 - i.e. about a 
factor 2 too high - constitutes therefore a severe anomaly 
with respect to the usual classification adopted for 0-type 
stars'^. In addition, a very massive star is also expected 
to be very luminous. From the observed V magnitude 
of HD 15558 and considering a reddening law with Ry 
= 3.1 and a color excess E(B - V) = 0.71'*, we derive 
an extinction Ay equal to 2.20, leading to an apparent 
deredde ned V magnitude equal to 5.70. Using the relation 
given by iHowarth fc Prinial l)l989l) , we obtain a bolomet- 
ric correction of -3.86, yielding a dereddened apparent 
bolometric magnitude equal to 1.84. Considerin g a dis- 
tance to IC1805 of 2.3 kpc l)Massev et al.lll995|) . we de- 
rive a log(Lf,oi/L0) equal to about 5.9. If we consider that 
the primary contributes about 5/6 of the total bolomet- 
ric luminosity, log(Lf,oi/Lo) reduces to about 5.82. This 
value is much too low to be compatible with a 150 
primary. Indeed, a confrontation of the bolometric lumi- 
nosity t o theoretical evolutionary tracks for vario us metal- 
licities (jSchaller et alJll993lSchaerer et al suggests 
masses not higher than about 70 Mq . If the distance we 
assumed here is correct, such a massive primary should be 
much more luminous than observed. 

On the other hand, an initially extremely massive 
and very hot early-type (i.e. spectral type 02) main 
sequence star will cool down during its evolution towards 
the giant luminosity cla ss. Fo r instance, according to the 
evolutionary models of ISchaller et al 

J ill 99 j) . a star of 
initial mass 120 M© starts its evolution with an effective 
temperature of 53000 K. When the effectiv e temperature 
reach es 38000 K, typical for an 05.5 giant l|Martins et al.l 
l2005l) . the stellar mass has decreased to about 85 Mq. 
Moreover, at this evolutionary stage, the hydrogen 
surface abundance should already be reduced, whilst 
helium should be enhanced. The main conclusion here is 
that a very massive star would have to lose a substantial 

^ We rely on the fact that the confidence interval derived 
by our analysis shows that there is a 65 % probability that 
the minimum mass of the primary lies between about 100 and 
200 Mq . This does not rule completely out a situation where 
this mass is lower than 100 Mq . 

^ We note also that a mass of 50 Mq for the 07V secondary is 
somewhat excessive. However, considering the error bar on the 
mass of the secondary, the 1-a confidence i nterval is still com- 
patibl e with the typical values proposed bv lHowarth fc Prinial 
il989l) for an 07V star, i.e. 36 Mq. Consequently, we consider 
that the mass of the primary is the main issue worth addressing 
here. 

* We obtained the co lor excess E(B - V) u sing the observed 
(B - V) =0.5 given bv lMassev et alJ 1^9^ and the intrinsic 
color (B - V)o based on the spectral type. 



fraction of its initial mass before reaching the effective 
temperature of an 05.5 giant. This would then imply an 
extremely large initial mass for the primary of HD 15558. 

Let us consider an alternative scenario, where 
HD 15558 is not a binary but a hierarchical triple sys- 
tem. The primary may be constituted of a yet unrevealed 
close binary system. In this scenario, as the mass of the 
primary is estimated on the basis of the motion of the 
secondary, the primary object - i.e. the hypothetical close 
binary - would appear to be a massive object whose mass 
is the sum of the masses of two stars. Provided the spectral 
type derived from our spectra for the primary is typical of 
the two stars constituting the close binary, we are possibly 
observing the composite spectrum of two similar 0-stars, 
in addition to that of the secondary object whose spectral 
type should be 07. This scenario offers the possibility to 
explain the unexpected high mass of the primary object, 
and to reconcile its mass with its spectral type. However, 
the triple system scenario still needs to solve the following 
issue: our spectral time series did not reveal any binary 
motion on a time-scale of a few days. If the primary is in- 
deed constituted of two stars, the short period orbit might 
be seen under a very low inclination angle. We might also 
consider the possibility that the time-scale of this orbit is 
a few weeks or so, and therefore poorly sampled by our 
spectral time series. However, we note that the radial ve- 
locity curve plotted in Fig.^does not present any strong 
dispersion likely to be due to an orbital motion on a time- 
scale significantly shorter than the main period of 442 d. 
In addition, we might expect some additional X-rays to be 
produced by an interaction between the winds of the two 
stars constituting the primary. High quality data obtained 
for instance with XMM-Newton are strongly needed to in- 
vestigate the X-ray emission from HD 15558 and discuss 
its origin in detail. In summary, even though our data do 
not provide any strong evidence supporting the triple sce- 
nario, this scenario may reconcile the mass derived for the 
primary and its spectral type. 

5.3. An open duster with unusually large stellar 
masses? 

The existence of such a massive star is worth discussing 
in the context of the massive star population of IC 1805, 
which should not be addressed without considering the 
presence of the most evolved object known in the vicin- 
ity: the microquasar LSI -|-61°303. This high mass X-ray 
binary producing a collimated relativistic jet consists of 
a Be star and of a compact object whose nature (neu- 
tron star or black hole) is not yet established l|Massi 
|2004|) . LSI-h61°303 is beheved to have been ejected out 
of IC 1805 during the supernova explosion of the initially 
most massive component of the binary l|Mirabel et alJ 
I2OO4I) . If the progenitor of LSI -f 61°303 was formed at the 
same epoch as the other 0-stars in IC 1805, among which 
is HD 15558, its primary component may have been more 
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massive than the primary of HD 15558. Alternatively, if 
LSI +61°303 was part of an older population of stars, its 
supernova explosion might have triggered the formation 
of the current population of 0-type stars. 

Two additional stars of this cluster are believed to have 
very large stellar masses, namely HD 15570 and HD 15629 
(see^Hcrrcro ct al. 2000), even though their masses have 
up to now only been estimated through model atmosphere 
fits. Unfortunately, as these two stars are probably single, 
an indepcndant mass determination through the study of 
a binary motion is unlikely. 

This series of presumably very massive objects sug- 
gests that IC 1805 harbours a population of particu- 
larly massiv e stars as compared to other open clusters. 
According to lMassev et al.l ()l995(l . a large number of very 
massive stars in an open cluster may be explained by its 
youth. Indeed, the age of the massive star population in 
IC 1805 was estimated to be 2 ± 1 Myr. Its most massive 
members have therefore not yet evolved into compact ob- 
jects. 

5.4. Multiplicity in IC1805 

The multiplicity of mas sive stars in young open clusters is 
a crucial question. The ICarcia jT^ ermilliod ( 2001) pro- 
posal of a binary frequency of 80 % for 0-stars suggests 
that the star formation process that was at work in IC 1805 
favored the formation of massive binary systems. However, 
in our s pectroscopic study of three mid-O-type stars of this 
cluster l)R,auw fc De Becked Eool . we showed that only 
one system is a spectroscopic binary (BD -1-60° 497) , whilst 
the other two are most probably single stars (BD -1-60° 501 
and BD -1-60° 513), leading t o the concl usion that the bi- 
nary frequency claimed by ICarcfa iTl vIcrmilliod ( 2001) 
might be overestimated. 

In Table|3 we summarize our present knowledge of 
the multiplicity of the 0-star members of IC 1805. At 
this stage, spectroscopic monitoring confirmed the bina- 
rity of only two members, i.e. HD 15558 (this paper) and 
BD +60° 497 (Paper I; Hillwig ct al. 200^, and our inves- 
tigations revealed no indication of binarity for four of them 
(BD-h60°501, BD-H60°513, HD 15570 and HD 15629). 
For the four remaining 0-type members, the multiplicity 
remains an open question. Consequently, at this stage all 
we can say is that the 0-star binary frequency in IC 1805 
should be of at least 20 % but is most probably not more 
than 60%. 

6. Summary and conclusions 

We have presented the results of an intensive spectroscopic 
study of the brightest massive stars in the young open 
cluster IC 1805: HD 15570, HD 15629 and HD 15558. For 
HD 15570 and HD 15629, the RVs do not present any sig- 
nificant trend attributable to binary motion on a time 
scale of a few days, nor from one yea r to the next. This is 
in line with the results obtained bv iHillwig et all ll2006t> 
who performed a search for RV variation on time-scales 



Table 3. 0-type star content of IC 1805. The first colurn n 
gives the star number following IVasilevskis et all l) 19651) . 
For the multiplicity, 's' means that our investigations did 
not reveal any indication of binarity, and '?' a lack of spec- 
troscopi c monitoring. The ref erences for the spe ctral types 
arej^l) iMassev et a1.l iTmfli. (2) Paper I, (3) lUnderhilJ 
l|l967l) . (4) llshidal ill 97(1) . (5) this study (Notes: * the 
spectral type of the four stars that are not studied in 
Paper I or in this paper may be uncertain as both O- and 
B-types were proposed by various authors, ** the spec- 
tral type of the secondary should range between 08.5 and 
09.5, *** IUnderhilil l|l967l) reported double hues for this 
star.) 



# 


ID 


Sp. Type* 


Ref. 


Status 


21 




09.5V((f)) 


1 


? 


104 


BD +60° 497 


06.5V((f)) 


2 


SB2 






+ 09V((f))** 






112 


BD +60° 498 


09V 


3 


?*•• 


113 




09Ve 


4 


? 


118 


BD +60° 499 


09.5V((f)) 


1 


? 


138 


BD +60° 501 


07V((f)) 


2 


s 


148 


HD 15558 


05.5IH(f) 


5 


SB2 






+ 07V 






160 


HD 15570 


04If+ 


5 


s 


192 


HD 15629 


05((f)) 


5 


s 


232 


BD+60° 513 


07.5V((f)) 


2 


s 



of a few days for some 0-type stars of IC 1805 including 
HD 15629 and HD 15570. 

In the case of HD 15558, we have derived an SBl 
orbital solution with significantly refined parameters 
as com pared to those obtained by ICarmanv &: Massevl 
The system appears to be eccentric (e~0.4) and 
we obtain a period of 442 ± 12 d. A careful inspection of 
the spectra obtained close to the extrema of the radial 
velocity curve reveals the presence of the companion in 
the profiles of the Hei A 4471, Hen A 4542, Civ A 5812 
and He I A 5876 lines. We have simultaneously fitted 
Gaussians to the profiles of the Hei A 4471 and Heii A 
4542 lines in order to separate the primary and secondary 
components. The determination of the equivalent width 
of these two lines allowed us to derive 05.5 and 07 
spectral types respectively for both stars. Considering in 
addition that the He ii A 4686 line of the primary has a 
P-Cygni profile whilst the same line is in absorption in 
the case of the secondary, along with the fact that we 
do not clearly observe the Niii AA 4634-4641 lines in 
emission for the secondary, we propose that HD 15558 is 
an 05.5111(f) + 07V binary. 

We estimated the radial velocities of the secondary of 
HD 15558 following two approaches: (1) a simultaneous fit 
of line profiles and (2) a disentangling method. Both tech- 
niques allowed us to determine minimum masses of the 
order of 150 ± 50 and 50 ± 15 respectively for the pri- 
mary and the secondary. We also obtained the first SB2 



12 



M. De Becker et al.: Early-type stars in the young open cluster IC 1805 



orbital solution for HD 15558. Although we note that the 

quality of this SB2 solution is rather poor, our data point 
clearly to a rather high mass ratio (about 3), leading to 
an extreme minimum mass for the primary. Our results 
require independent validation using an improved disen- 
tangling procedure. 

The main problem in considering our results is to rec- 
oncile the very extreme mass of the primary with its spec- 
tral type. It is indeed unlikely that a very massive main- 
sequence star could cool down enough during its evolution 
to become an 05.5 giant. A possible scenario can however 
be considered where HD 15558 is not a binary but a triple 
system. The primary may be a yet unrevealed close bi- 
nary system. In this case, as the mass of the primary is 
estimated on the basis of the motion of the secondary, 
the primary object - i.e. the hypothetical close binary - 
would appear to be a massive object whose mass is the 
sum of the masses of two stars. Even though our data do 
not provide any evidence for this scenario, we estimate 
that at this stage it should not be rejected, and that it 
could constitute a valuable working hypothesis for future 
investigations concerning this system. 

From our new and previously published results, we 
briefly address the question of the multiplicity of the early- 
type stars in IC 1805. The binary frequency among 0-stars 
should be of at least 20%, since out of 10 0-stars only 2 
are confirmed binaries, and should not exceed 60%. We 
therefore conclude that the previously claimed binary fre- 
quency of 80 % was overestimated. 
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Table A.l. Observing runs used for the line profile vari- 
ability study of HD 15570, HD 15629 and HD 15558. The 
first and second columns give the name of the campaign 
as used in the text as well as the instrumentation used. 
The next columns are the number of spectra obtained, the 
time elapsed between the first and the last spectrum of the 
run, the natural width of a peak of the power spectrum 
taken as 1/AT, and the mean signal-to-noise ratio of each 
data set. 



Obs. run 


Telescope 


N 


AT 




S/N 








(d) 


(d-) 




HD 15570 


Sept. 2000 


OHP/1.52m 


7 


9.99 


0.10 


390 


Sept. 2001 


OHP/1.52m 


5 


7.02 


0.14 


360 


Sept. 2002 


OHP/1.52m 


10 


14.94 


0.07 


300 


Oct. 2003 


OHP/1.52m 


7 


17.97 


0.06 


300 


Oct. 2004 


OHP/1.52m 


7 


9.90 


0.10 


260 


Oct. 2004 


SPM/2.10m 


4 


2.07 


0.48 


180 


HD 15629 


Sept. 2002 


OHP/1.52m 


9 


12.99 


0.08 


300 


Oct. 2003 


OHP/1.52m 


4 


17.98 


0.06 


280 


Oct. 2004 


OHP/1.52m 


7 


9.90 


0.10 


260 


HD 15558 


Sept. 2000 


OHP/1.52m 


12 


10.97 


0.09 


530 


Sept. 2001 


OHP/1.52m 


12 


7.02 


0.14 


440 


Sept. 2002 


OHP/1.52m 


11 


14.92 


0.07 


430 


Oct. 2003 


OHP/1.52m 


7 


17.95 


0.06 


360 


Oct. 2004 


OHP/1.52rn 


5 


7.00 


0.14 


450 



Appendix A: Observations and data reduction 

Spectroscopic observations were collected at the 
Observatoire de Haute-Provence (OHP, France) during 
several observing runs from 2000 to 2004 for HD 15570 
and HD 15558, and from 2002 to 2004 for HD 15629. AU 
spectra were obtained wit h the Aurelie spect rograph fed 
by the 1.52 m telescope (jCillet et all Il994l) . using the 
same setup as described in Paper I. Our collection of 
spectra is described in Table lA.l1 

We also obtained several spectra of HD 15570 with 
the echelle spectrograph mounted on the 2.1m telescope 
at the Observatorio Astronomico Nacional of San Pedro 
Martir (SPM) in Mexico, with exposure times ranging 
from 10 to 20 minutes. The instrument covers the spectral 
domain between about 3800 and 6800 A . The detector 
was a Site CCD with 1024 x 1024 pixels of 24 /xm^. The 
slit width was set to 200 /im corresponding to 2 arcsec on 
the sky. The data were reduced using the echelle package 
available within the MIDAS software. After adding 
consecutive spectra of a given night to reach higher 
signal-to-noise ratios, at the expense of time resolution, 
we obtained 4 spectra of HD 15570. 



We obtained 17 additional s pectra of HD 15558 with 
the Elodie echelle spectrograph ijBaranne et al ] IT99i fed 
by the 1.93 m telescope at the Observatoire de Haute 
Provence between March 2003 and February 2005 to mon- 
itor a complete orbital period. The exposure time of 
each of these spectra was 90 minutes. This spectrograph 
uses a combination of a prism and a grism as a cross- 
disperser, with a blaze angle of 76°. The resolving power 
achieved is about 42000 between 3906 and 6811 A in a 
single exposure, and the detector is a Tkl024 CCD with 
24 /im X 24 /im pixels. The Elodie data consist of single 
spectra distributed over 67 orders. Due to pointing con- 
straints specific to the 1.93 m telescope, no echelle spec- 
tra were obtained between April 2004 and July 2004. We 
filled that gap with 6 observations with the Aurelie spec- 
trograph mounted on the 1.52 m telescope. Two spectra 
out of the 6 were obtained using the same grating as 
described in the previous paragraph (between 4455 and 
4890 A ) . The four remaining spectra were obtained using 
a 1200 1/mm grating providing a resolving power of about 
16000 in the blue range, with a reciproqual dispersion of 8 
Anim^^ (between 4455 and 4680 A). The exposure time 
of these 6 Aurelie spectra was 60 minutes. 

For Aurelie data, we adopted the same reduction pro- 
cedure as in Paper I. Elodie data were reduced using the 
standard on-line automatic treatment implemented at the 
OHP. All spectra were normalized using splines calculated 
on the basis of properly chosen continuum windows. 

Appendix B: Radial velocity measurements 

The mean RVs of HD 15570 and HD 15629 obtained for the 
various observing runs are collected in Table lBJl In the 
case of HD 15558, all the RVs used to obtain the orbital 
solution are listed in Table lR2l All the RVs considered 
in this pap er were estimated u sing the rest wavelengths 
provided bv lConti et alJ l)l977j) . 
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Table A. 2. Description of the data of HD 15558 obtained 
during the long-term monitoring using the 1.93 and the 
1.52 m telescopes. The first column gives the date of the 
observation. The instrumentation used to obtain the spec- 
trum is provided in the second column. The next columns 
give the resolving power of the instrumentation used, and 
the signal-to-noise ratio. 



Date 


Telescope /Instr. 


A/AA 


S/N 


2003/09/04 


1.93 m/Elodie 


42000 


170 


2003/10/04 


1.93m/Elodie 


42000 


160 


2003/10/21 


1.93 m/Elodie 


42000 


100 


2004/01/07 


1.93 m/Elodie 


42000 


180 


2004/01/15 


1.93 m/Elodie 


42000 


100 


2004/02/14 


1.93 m/Elodie 


42000 


130 


2004/03/11 


1.93 m/Elodie 


42000 


110 


2004/04/04 


1.93 m/Elodie 


42000 


110 


2004/04/10 


1.52 m/Aurelie 


8000 


330 


2004/05/06 


1.52 m/Aurelie 


16000 


160 


2004/06/09 


1.52 m/Aurelie 


16000 


290 


2004/06/28 


1.52 m/Aurelie 


16000 


360 


2004/07/15 


1.52 m/Aurelie 


8000 


550 


2004/07/26 


1.52 m/Aurelie 


16000 


360 


2004/08/15 


1.93 m/Elodie 


42000 


100 


2004/08/26 


1.93 m/Elodie 


42000 


120 


2004/09/06 


1.93 m/Elodie 


42000 


140 


2004/10/07 


1.93 m/Elodie 


42000 


130 


2004/11/09 


1.93 m/Elodie 


42000 


60 


2004/11/17 


1.93 m/Elodie 


42000 


140 


2004/12/03 


1.93 m/Elodie 


42000 


140 


2005/01/27 


1.93 m/Elodie 


42000 


100 


2005/02/16 


1.93 m/Elodie 


42000 


70 



Table B.l. Radial velocity of the main He and N lines 
for both stars measured on OHP spectra (expressed in 
kms^^). The mean radial velocity along with its 1-a stan- 
dard deviations are provided for each observing run. All 
individual RVs can be found in the WEBDA data base at 
|http : //obswww . unige . ch7webda| 



HD 15570 



Data set 


He I 


Hon 


Nlll 


Nni 




A 4471 


A 4542 


A 4634 


A 4641 


Sept. 2000 


-69.1 ± 8.2 


-46.7 ± 8.7 


-56.1 ± 3.9 


-67.0 ± 2.8 


Sept. 2001 


-65.1 ± 10.6 


-51.1 ± 3.2 


-58.6 ± 3.5 


-70.3 ± 2.0 


Sept. 2002 


-60.5 ± 13.2 


-50.0 ± 5.0 


-56.8 ± 4.6 


-70.0 ± 3.3 


Oct. 2003 


-59.3 ± 15.7 


-46.2 ± 3.9 


-55.6 ± 4.5 


-67.6 ± 5.4 


Oct. 2004 


-72.7 ± 8.3 


-47.8 ± 7.0 


-58.0 ± 4.7 


-70.6 ± 3.7 


All 


-65.0 ± 12.3 


-48.3 ± 5.9 


-56.9 ± 4.2 


-69.1 ± 3.8 






HD 15629 




Data set 


He I 


Hon 


Nlll 


Niii 




A 4471 


A 4542 


A 4634 


A 4641 


Sept. 2000 










Sept. 2001 










Sept. 2002 


-60.6 ± 3.6 


-46.9 ± 4.5 


-59.9 ± 8.8 


-72.1 ± 6.8 


Oct. 2003 


-55.3 ± 4.6 


-41.0 ± 1.6 


-54.6 ± 3.2 


-71.6 ± 2.6 


Oct. 2004 


-58.2 ± 5.9 


-50.1 ± 4.6 


-63.5 ± 8.3 


-74.6 ± 6.5 


All 


-58.7 ± 4.9 


-46.8 ± 5.2 


-60.1 ± 8.2 


-72.9 ± 6.0 



Table B.2. Radial velocities obtained from our time se- 
ries of HD 15558. The second and third columns give the 
heliocentric Julian day (-2 450 000) and the orbital phase 
following the parameters provided in Tabled (left part). 
The next columns provide the radial velocities obtained 
for lines that were selected for the determination of the 
SBl orbital solution. The column labelled 'Mean' contains 
the mean of the radial velocities obtained for the He ii and 
the Niii lines quoted in the fourth and fifth columns. All 
RVs are expressed in kms~^. The last column gives the 
weight (W) attributed to our measurements to calculate 
the orbital parameters, depending on the spectral resolu- 
tion of the instrument and the signal-to-noise ratio of the 
spectra. 



# H,]D (j) Hen Niii Mean W 

A 4542 AA 4634, 



4641 


1 


1810. 


,640 


0, 


.034 


-85.9 


-96.1 


-91.0 


0.25 


2 


1810. 


,658 


0, 


.034 


-88.0 


-117.0 


-102.5 


0.25 


3 


1811. 


,619 


0, 


.037 


-86.7 


-109.0 


-97.9 


0.25 


4 


1811. 


,665 


0, 


.038 


-88.8 


-101.5 


-95.1 


0.25 


5 


1812. 


,660 


0, 


.039 


-87.5 


-112.6 


-100.0 


0.25 


6 


1813. 


,651 


0, 


.041 


-87.8 


-118.8 


-103.3 


0.25 


7 


1814. 


,643 


0, 


.043 


-91.6 


-118.9 


-105.2 


0.25 


8 


1815. 


,645 


0, 


.046 


-88.9 


-110.2 


-99.6 


0.25 


9 


1818. 


,641 


0, 


.053 


-86.4 


-110.6 


-98.5 


0.25 


10 


1819. 


,593 


0, 


.055 


-86.8 


-109.9 


-98.3 


0.25 


11 


1820. 


,648 


0, 


.057 


-83.0 


-118.1 


-100.6 


0.25 


12 


1821. 


,609 


0, 


.059 


-91.2 


-111.2 


-101.2 


0.25 


13 


2163. 


,580 


0, 


.832 


-1.2 


-12.1 


-6.7 


0.25 


14 


2164. 


,597 


0, 


.834 


-4.3 


-10.4 


-7.3 


0.25 


15 


2164. 


,660 


0, 


.835 


-0.9 


-18.8 


-9.9 


0.25 


16 


2165. 


,577 


0, 


.837 


-3.5 


-23.7 


-13.6 


0.25 


17 


2165. 


,590 


0, 


.837 


-0.1 


-12.9 


-6.5 


0.25 


18 


2167. 


,550 


0, 


.841 


-1.9 


-7.4 


-4.6 


0.25 


19 


2167. 


,564 


0, 


.841 


-6.7 


-14.8 


-10.7 


0.25 


20 


2168. 


,601 


0, 


.843 


-1.8 


-17.3 


-9.5 


0.25 


21 


2168. 


,642 


0, 


.844 


-6.0 


-18.9 


-12.4 


0.25 


22 


2169. 


,581 


0, 


.846 


-4.2 


-16.8 


-10.5 


0.25 


23 


2170. 


,580 


0, 


.848 


-7.2 


-22.1 


-14.6 


0.25 


24 


2170. 


,603 


0, 


.848 


-8.1 


-23.9 


-16.0 


0.25 


25 


2518. 


,617 


0, 


.634 


-27.5 


-28.5 


-28.0 


0.25 


26 


2520. 


,593 


0, 


.639 


-14.8 


-20.2 


-17.5 


0.25 


27 


2523. 


,580 


0, 


.646 


-11.5 


-32.8 


-22.2 


0.25 


28 


2524. 


,510 


0, 


.648 


-13.8 


-28.0 


-20.9 


0.25 


29 


2527. 


,531 


0, 


.655 


-12.4 


-21.6 


-17.0 


0.25 


30 


2528. 


,499 


0, 


.657 


-14.5 


-26.7 


-20.6 


0.25 


31 


2529. 


,529 


0, 


.659 


-10.5 


-21.6 


-16.1 


0.25 


32 


2531. 


,503 


0, 


.664 


-4.9 


-21.8 


-13.4 


0.25 


33 


2532. 


,516 


0, 


.666 


-8.3 


-20.6 


-14.5 


0.25 


34 


2532. 


,644 


0, 


.666 


-7.9 


-27.2 


-17.5 


0.25 


35 


2533. 


,532 


0, 


.668 


-12.6 


-24.2 


-18.4 


0.25 


36 


2886. 


,627 


0, 


.466 


-37.0 


-51.5 


-44.2 


1.00 


37 


2916. 


,592 


0, 


.534 


-29.4 


-62.0 


-45.7 


1.00 


38 


2916. 


,558 


0, 


.534 


-24.7 


-39.1 


-31.9 


0.25 


39 


2918. 


,570 


0, 


.538 


-25.7 


-52.6 


-39.2 


0.25 


40 


2919. 


,561 


0, 


.541 


-26.6 


-47.8 


-37.2 


0.25 


41 


2922. 


,566 


0, 


.547 


-24.4 


-34.4 


-29.4 


0.25 


42 


2926. 


,573 


0, 


.556 


-19.6 


-34.8 


-27.2 


0.25 


43 


2928. 


,593 


0, 


.561 


-24.7 


-39.4 


-32.0 


0.25 
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# HJD 4, He II Niii Mean W 

A 4542 AA4634, 
4641 



45 


2934.481 


0.574 


-13.7 


-47.1 


-30.4 


1, 


.00 


46 


3012.476 


0.751 


-14.0 


-31.5 


-22.7 


1, 


.00 


47 


3020.329 


0.768 


-19.2 


-20.7 


-20.0 


1, 


.00 


48 


3050.375 


0.836 


-11.2 


24.6 


-17.9 


1 


.00 


49 


3076.317 


0.895 


-16.7 


-34.3 


-25.5 


1 


.00 


50 


3100.318 


0.949 


-43.1 


-62.1 


-52.6 


1 


.00 


51 


3106.340 


0.963 


-28.8 


-47.1 


-38.0 





.25 


52 


3132.374 


0.022 


-63.5 


-85.2 


-74.3 





.50 


53 


3165.550 


0.097 


-88.2 


114.4 


-101.3 





.50 


54 


3184.590 


0.140 


-82.7 


104.4 


-93.6 





.50 


55 


3201.593 


0.178 


-72.1 


98.7 


-85.4 


0, 


.25 


56 


3212.544 


0.203 


-65.3 


-94.9 


-80.1 


0, 


.50 


57 


3232.579 


0.248 


-67.1 


-87.5 


-77.3 


1, 


.00 


58 


3244.384 


0.275 


63.0 


-93.0 


-78.0 


1 


.00 


59 


3254.630 


0.298 


-46.3 


-78.0 


-62.2 


1, 


.00 


60 


3286.477 


0.370 


-43.8 


-60.1 


-52.0 


1, 


.00 


61 


3289.570 


0.377 


-46.3 


-62.9 


-54.6 


0, 


.25 


62 


3290.553 


0.379 


-35.6 


-53.3 


-44.4 





.25 


63 


3294.603 


0.388 


-32.3 


46.6 


39.4 





.25 


64 


3295.584 


0.390 


-37.6 


-44.2 


-40.9 


0, 


.25 


65 


3296.569 


0.393 


-44.0 


-56.0 


-50.0 





.25 


66 


3318.587 


0.442 


-28.5 


-41.8 


-35.1 





.10 


67 


3326.508 


0.460 


33.4 


-56.5 


45.0 


1 


.00 


68 


3343.487 


0.499 


-27.9 


-54.4 


41.2 


1 


.00 


69 


3398.308 


0.623 


-9.4 


-33.8 


-21.6 


1, 


.00 


70 


3418.295 


0.668 


-3.6 


-24.5 


-14.0 


0, 


.10 



